Anisotropies in the Motions and Positions of the Galactic 

Globular Clusters 



F.D.A. Hartwick 
Department of Physics and Astronomy, 
University of Victoria, Victoria, BC, Canada, V8W 3P6 

ABSTRACT 

The velocity elhpsoid for 38 globular clusters with [Fe/H]< —1.0 is derived 
and shown to be significantly anisotropic with major axis directed towards low 
Galactic latitude. Principal axes of the spatial distribution of different groups of 
clusters are derived and compared with the velocity ellipsoid. The metal poor 
cluster spatial distribution is significantly flattened along an axis which coincides 
within the uncertainties with the major axis of the velocity ellipsoid. Given the 
observed steep age-metallicity relation for metal poor clusters, one speculative 
interpretation of the data is that an initially flattened filament underwent a 
relatively rapid initial transverse collapse forming satellite galaxies and metal 
poor globular clusters while the protogalaxy collapsed and assembled more 
slowly along the filament acquiring and/or redistributing angular momentum in 
the process. 



Subject headings: Galaxy: halo — Galaxy: formation — globular clusters: 
general 



1. Introduction 



Globular clusters provide a unique probe of the earliest phases of the evolution of 
the Galaxy. From the early work of Kinman (1959), Zinn (1985), and others it is clear 
that the clusters can be divided into at least two groups - the spatially extended, metal 
poor blue clusters and the more centrally concentrated towards the Galactic center, metal 
rich red clusters. This work is concerned primarily with the former group, and the main 
motivation is the recent availability of space motions for a significant number of these 
clusters as published by Dinescu et al. (1999). Dinescu (private communication) has also 
provided this author with data for 4 additional clusters. These data are used to determine 
the velocity ellipsoid for the metal poor group in order to allow a comparison with the 



- 2 - 



spatial distribution. In a previous paper (Hartwick 2000, hereafter H2000) the spatial 
distribution of various halo samples was examined. There it was found that the metal poor 
group of globular clusters formed an oblate distribution whose minor axis was directed to 
a low Galactic latitude. This was in contrast to the metal rich group which exhibited a 
triaxial distribution with minor axis toward the Galactic pole. We re-examine the spatial 
distributions with modified samples of clusters and compare the results with the kinematics. 



2. The Velocity Dispersion Tensor of the Halo Clusters 

Danescu et al. (1999) have compiled space motions for 38 globular clusters, and 35 of 
these have [Fe/H] values less than —1.0. Since that time 3 more halo clusters (NGC6522, 
NGC7006, & Pal 13) have been added to the compilation. This sample of 38 clusters is 
the basis for the following kinematical investigation. The tabulated UVW components of 
space motion (corrected solar pecuhar motion and for Galactic rotation of 220 km sec~^) 
are referred to to the Galactic center where T4 is the radial component in the Galactic 
plane directed positive away from the center, is the tangential component positive in 
the direction of Galactic rotation and is the vertical component directed positive in the 
direction of the north Galactic pole. 

The average of each of these components for the entire sample was removed before 
forming the following six moments: ^0^0, K-^^, K-Vz, and V^Vz. These six 

moments are the components of the symmetrical velocity dispersion tensor. The eigenvalues 
of the matrix were then determined in the standard way and are given below. The units 
of the eigenvalues denoted by a arc km sec~^ and the one sigma uncertainties throughout 
were computed by the bootstrap method. 

Clusters with [Fe/H]< -1.00 (n=38) 
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As a check on the sensitivity of the above result to the [Fe/H] cutoff, we repeated the 
procedure for the 26 clusters with [Fe/H]< —1.5 and find: 



Clusters with [Fe/H]< -1.50 (n=26) 
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a, = 122l«f, l^mtlli. b^-S2.SXl%l (2) 
(Jc = 94.71^7^^ I = 60.8lii:i b = -3.481^?:^ 
n = 49.311:1 



In both of the above solutions the mean streaming motions in the vr and z directions 
differed from zero by less than 1.5 sigma. Also the above values of are independent of 
distance from the Galactic center. 

Asymmetric error ellipsoids in the components of space motion can be expected due 
to relatively small errors in radial velocity and relatively large errors in proper motion. 
In order to ensure that such effects were not unduly influencing solution (1) the errors 
in radial velocity, proper motion and distance given by Dinescu et al. (1999) were used 
as dispersions in assumed gaussian distributions and solutions were performed on space 
motions determined by the propagation of randomly chosen uncertainties in each of the 
above quantities. The mean and standard deviation of the major axis from 100 trials with 
38 clusters is (7„ = 167 ± 7.98 towards Z = 332 ± 6.28 and b = -5.94 ± 7.64. Results for the 
other two axes are ai, — 122 ± 10.9 and ctc = 108 ± 10.1. This independent assessment of 
uncertainties thus strengthens the case for the generally prolate velocity ellipsoid given by 
solution (1). 

An examination of the spatial distribution of the clusters from solution (1) shows that 
only 7 of the 38 clusters were in the hemisphere beyond the Galactic center. In order to 
assess the effect of this spatial imbalance on our result, kinematic solutions were made using 
the above 7 clusters combined with 7 independently chosen clusters from the hemisphere 
containing the sun. Thirty solutions were made with 14 clusters each and the median 
results for each aspect of the major axis are: aa = IGGlg^, / = 348^7, and b = —13^1^. 
While the one sigma uncertainties of solution (1) do overlap with the above result, perhaps 
the most important implication of this exercise is a clear need for more cluster space motion 
determinations. 



3. The Spatial Distribution of the Globular Clusters 

The spatial distribution of the outlying satellites and the globular clusters was 
determined in H2000. There it was found that the 99 metal poor globular clusters 
([Fe/H]< —1.0) formed an oblate distribution with minor axis directed towards 
I — 307°, 6 = —3.3° quite different from the short axis of the distribution of the satellite 
galaxies. The result from the kinematics described above prompted a re-examination of this 
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result. It is believed that the pole of the orbit of the Sagittarius system is highly inchned to 
that of the remaining outlying satellite galaxies (Lynden-Bcll & Lyndcn-Bcll 1995). As well 
these authors tentatively assign the 6 clusters M54, NGC2419, Arp 2, Pal 2, Terzan 7 and 
Terzan 8 to the Sagittarius stream. Given that the present orbit of Sagittarius is probably 
not the one that it had at formation (c.f. Zhao 1998, also Johnston et al. 2002) we removed 
the above clusters from the original sample and performed new solutions. (Note that the 
Harris (1996) web-based compilation contains 147 clusters of which 99 have [Fe/H]< —1.0. 
CuUing 5 clusters from this sub-group still leaves us with 95% of the original sample. 
Furthermore, none of these five clusters was included in the above kinematic solutions.) 
The analysis procedure followed here is shghtly different from that in H2000. Here the 
tensor components XiXi, Xiyl etc were calculated both without weighting and with weighting 
by l/d,j where the Xi, iji, and Zi are the projections of dj, the distance from the Galactic 
center (Rq = 8 kpc assumed). Both the weighted and unweighted solutions arc given since 
the true distribution probably lies somewhere in between. The eigenvalues of the resulting 
matrix are denoted Ca, e;,, and Cc and have units of kpc. 



Cluster Sample with [Fe/H]< -1.0 (weighted) (n=94) 
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Cluster Sample with [Fe/H]< -1.0 (unweighted) (n=94) 
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As a check on the robustness of the result, solutions were also made for the 60 remaining 
clusters with [Fe/H]< -1.5. 
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Cluster Sample with [Fe/H]< —1.5 (unweighted) (n=60) 
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We note that for both samples the short axis of the (nearly oblate) cluster distribution is 
very similar to the long axis of the kinematic solution. In order to insure that our result 
is not affected by incompleteness or uncertain absorption corrections, solutions were also 
performed for only those clusters whose absolute Galactic latitudes were larger than 10°. 



Cluster Sample with [Fe/H]< -1.0 and |6| > 10 (weighted) (n=66) 
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Cluster Sample with [Fe/H]< —1.0 and \b\ > 10 (unweighted) (n=66) 
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The quoted uncertainties in the above solutions were determined by the bootstrap 
method. It is instructive to consider the effects on the solution of errors in the individual 
distances to each cluster. Randomly chosen errors drawn from a gaussian distribution with 
a sigma of 10% of each distance were propagated through and the mean and standard 
deviation of the minor axis of the spatial distribution from 100 trials with 66 clusters was 
found to be Cc = 4.67 ± 0.0949 towards Z = 327± 1.28 and b = -9.91 ± 1.02. The amplitudes 
of the other two axes are Ca — 7.98 ± 0.146 and = 7.11 ± 0.148. The good agreement 
with solution (7) strengthens the case for the generally oblate spatial distribution given by 
solution (7). 

For illustration, the distribution of clusters in this new coordinate system (solution (7)) 
is shown in Fig. 1 where x' are projections of Galactocentric distance along the major axis 
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and z' are projections along the minor axis. As summarized in Table 1, the short axes of 
the spatial distributions are remarkably similar to the major axes of the kinematic sohitions 
of the previous section. Furthermore, the minor axis of the satellite galaxy distribution is 
also quite similar. 

In H2000 the spatial distribution of the 34 most metal rich globular clusters 
([Fe/H]> —0.7) was also determined. The distribution was found to be triaxial with the 
long axis directed slightly off (~ 17°) the sun-center line and the short axis now pointing to 
high Galactic latitude. We repeat the process below but with one less cluster (Terzan 7). 



Cluster Sample with [Fe/H]> -0.7 (weighted) (n=33) 
e„ = 2.02^°:^} I = 338tti:^ b = -6.09^^:1 

e, = 1.5lt',ill I = GS.ltirr b = OASOtlH (9) 

1-65.2 
-60.2 

c/a = 0.38t',Z 



Cluster Sample with [Fe/H]> -0.7 (unweighted) (n=33) 
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A triaxial distribution is still present, but now the major axis is ~ 26° off of the sun-center 
line while the minor axis is directed to high latitude. Again it seems reasonable to suggest 
that we are observing a globular cluster counterpart to the 'COBE' bar (c.f. Binney & 
Merrifield, 1998, p. 616). Should future work confirm the remarkably close alignment 
between the major axis of the spatial distribution and the major kinematic dispersion 
axis (see Table 1), an intimate connection between halo (metal poor clusters) and bulge 
(metal rich clusters) would be implied. One possible argument against such a connection 
is the recent suggestion that bars may be less common at high redshift than at z=0 (van 
den Bergh et al. 2002). Further elucidation will come from the determination of the space 
motions for these metal rich clusters. 



4. Discussion 



Both the velocity ellipsoid and the spatial distribution of the most metal poor Galactic 
globular clusters appear to be anisotropic. It is interesting that the major axis of the 
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velocity ellipsoid is directed towards low Galactic latitude, consistent with the direction of 
the short axis of the spatial distribution and not far from the pole of the OTitlying satellite 
galaxy distribution. A quantitative comparison for the different samples is given in Table 1 
where for reference we take the major axis of the kinematic solution (1). 

There would appear to be no straight forward explanation for the observations above. If 
the gravitational potential was spherical, one can argue that the direction of highest velocity 
dispersion would also be that of the longest spatial axis and vice versa. The above result 
then apparently rules out a spherical potential but would be consistent with one that was 
flattened in the direction of highest velocity dispersion. Recalling that the age-metallicity 
relation for the most metal poor globular clusters is extremely steep (vandenBerg 2000), 
a naive picture is that an initially flattened fllament underwent a rapid collapse tranverse 
to the long axis during which time satellite galaxies, the metal poor clusters, and the 
triaxial structure discussed in H2000 were formed. Simultaneously the Galaxy was being 
assembled by a slower collapse along the fllament during which time most of its angular 
momentum was being induced and/or redistributed. While triaxial structures arc formed 
in cold dark matter simulations, the scenario outlined above may be more consistent with a 
warm dark matter theory for structure formation where the lowest mass objects form from 
the fragmentation of caustics (c.f. Bode, Ostriker, & Turok 2001). 

The author wishes to thank Dr. Dana Dinescu for providing the latest compilation of 
cluster space motions in convenient form, and Ray Carlberg, Scott Tremaine, and Sidney 
van den Bergh for their helpful comments on an earlier draft. He also wishes to acknowledge 
flnancial support from an NSERC of Canada operating grant. 
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Table 1. Comparison of the Axes of the Various Samples 



Sample 
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samples. 



-9- 



vandenBerg, D. A. 2000, ApJS, 129, 315 

van den Bergh, S., Abraham, R.G.,Whyte, L.F., Merrifield, M.R., Eskridge, P.B., Progel, 
J.A. & Pogge, R. 2002, AJ, 123, 2913 

Zhao, H. 1998, ApJ, 500, L149 

Zinn, R. 1985, ApJ, 293, 424 



This preprint was prepared with the AAS IM^jX macros v4.0. 



-10- 




o 
o 



o 



o 

Q. 



o 

ID 



O 
O 



O 
O 



O 
ID 



O 



o 
o 



Fig. 1. — The spatial distribution of 66 globular clusters with [Fe/H]< —1.0 and \b\ > 10° 
in coordinates transformed to the major (x') and minor (z') axes of solution (7) illustrating 
the flattening of the distribution. For reference the minor axis is directed to I — 327°, 
b = -9.65°. 



